Annual change of Pareto index dynamically deduced from 
the law of detailed quasi-balance 



*Atushi Ishikawa 
Kanazawa Gakuin University, Kanazawa 920-1392, Japan 



Abstract 



Employing data on the assessed value of land in 1983 - 2005 Japan, we investigate the 
dynamical behavior in the high scale region of non-equilibrium systems. From the detailed 
quasi-balance and Gibrat's law, we derive a relation between the change of Pareto index and 
a symmetry in the detailed quasi-balance. The relation is confirmed in the empirical data 
nicely. 
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1 Introduction 

In physics, universal laws in nature are extracted by analyzing empirical data, and various 
phenomena are explained by the laws. In econophysics, economic phenomena are investigated 
by using this physical approach [1]. From this point of view, recently a major progress has been 
made with the research in distributions of personal income or company size. 

More than a hundred years ago, it is well known that a cumulative number N(> x) obeys a 
power-law for x which is larger than a certain threshold xq: 



This power-law and the exponent fj, are called Pareto's law and Pareto index, respectively [2] . 
Here x is wealth, income, profits, assets, sales, the number of employees and etc. 

Fujiwara et al. [5] have explained this phenomenon by using the law of detailed balance and 
Gibrat's law [6], which are observed in empirical data. The detailed balance is time-reversal 
symmetry: 




for x > Xq . 





(2) 
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1 Pareto's law is checked with high accuracy (See [3, 4] for instance). 



Here x\ and X2 are two successive incomes, profits, assets, sales, etc, and P\2(xi, X2) is a joint 
probability distribution function (pdf). Gibrat's law states that the conditional probability 
distribution of growth rate Q(R\x\) is independent of the initial value x\: 

Q{R\ Xl ) = Q(R) . (3) 

Here growth rate R is defined as the ratio R = X2/X1 and Q(R\x\) is defined by using the pdf 
P\(x\) and the joint pdf Plr(xi , R) as 

In the proof, they assume no model and only use these two underlying laws in empirical data. 

Below the threshold xq, it is also well known that the power-law is not observed [6, 7]. 2 In 
Ref. [11], by employing profits data of Japanese companies in 2002 and 2003 [12], Gibrat's law 
(3) is extended in the middle scale region as follows: 

Q(R\ Xl ) = d #-*+(*(>)-«+ in fori?>l, (5) 

Q(R\ Xl ) = dfl+MsoJ-a-inSi-i for R<1 (6) 

where 

u + ~ a_ ~ for x\ > xo , (7) 

u + ~ a_ 7^ for x m i n < xi < xo , (8) 

i+(x ) - i-(x ) = n ■ (9) 

By combining extended Gibrat's law (5) - (9) with the detailed balance (2), we have kinemati- 
cally derived the following distribution function in the high and middle scale region uniformly 

Pi(xi) = C Xl ~ {fl+l) e~ alv? % for x x > x min , (10) 

where a = (a+ + a_)/2. 

All these findings are enormously important for the progress of econophysics. Above deriva- 
tions are, however, valid only in the economic equilibrium where the detailed balance (2) holds. 
As the next step, the dynamics should be established. For this aim, we must investigate long- 
term economic data in which dynamical transitions are observed. Unfortunately, it is difficult 
to obtain personal income or company size data for a long period. 

In Ref. [13], Kaizoji has reported that the distribution of Japanese land prices has similar 
features with one of personal income and company size. The database of the assessed value of 
land is available on the Ministry of Land, Infrastructure and Transport Government of Japanese 
World-Wide Web site [14]. 

In this paper, we analyze the high scale region data on the assessed value of land in a 23-year 

period (1983 - 2005) in order to investigate the dynamical behavior. As a result, we obtain a 

dynamical equation which non-equilibrium economic systems should satisfy. 

2 It is pointed out that the reason is the breakdown of Gibrat's law [5]. The breakdown of Gibrat's law in 
empirical data is reported by Stanley's group [8]. Takayasu et al. [9] and Aoyama et al. [10] also report that 
Gibrat's law does not hold in the middle scale region by using data of Japanese companies. 
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2 Data of land prices and the law of detailed quasi-balance 



In Japan, land is a very important asset which is distinguished from building. The transition 
of land prices has a possibility to influence the economy. In fact, Japanese economy experienced 
bubble term (1986 - 1991) caused by the abnormal rise of land prices. 

The assessed value of land indicates the standard land prices evaluated by Ministry of Land, 
Infrastructure, and Transport. The investigation is undertaken on each piece of land assessed 
once a year. We employ the database of the assessed value of land covering the 23-year period 
from 1983 to 2005. The distributions of land prices are shown in Fig. 1, 2, 3, 4 and 5. The 
number of data points of land prices increases gradually, because the database only contains 
data points which exist in the 2005 evaluation. 

From Fig. 1-5, the power-law is confirmed in the high scale region. For each year, we 
estimate Pareto index ji in the range of land prices from 2 x 10 5 to 10 7 yen/m 2 where Pareto 
law holds approximately. Annual change of Pareto index \x from 1983 to 2005 is represented in 
Fig. 6. In this period, Pareto index has changed annually. This means that the system is not in 
equilibrium and the detailed balance (2) does not hold. Actually, the breakdown is observed in 
the scatter plot of all pieces of land assessed in the database (Fig. 7 and 8 for instance). There 
is no x\ <-> X2 symmetry in Fig. 7 and 8 obviously. On the other hand, the detailed balance 
(x\ <-> X2 symmetry) is observed approximately in Fig. 9 and 10 for instance. 

Here we make a simple assumption that the symmetry of the joint pdf Pi2(xi, X2) is repre- 
sented as a regression line fitted by least-square method as follows 



The detailed balance (2) has the special symmetry, 6 = a = 1. For each scatter plot, we measure 
6, a and the result is shown in Fig. 11. 

From this symmetry {ax\ 6 <-> X2), we extend the detailed balance (2) to 



In this paper, we call this law the detailed quasi-balance. 

3 Dynamically varying Pareto index 

In this section, we derive the relation between the change of Pareto index /x and in Eq. (11) 
by extending the detailed balance (2) to the detailed quasi-balance (12). We assume Gibrat's 
law in the high scale region, because the number of data points is insufficient to observe Gibrat's 
law. 

Due to the relation of P\2{x\,X2)dx\dx2 = P\r{x\, R)dx\dR under the change of variables 
from (xi,X2) to (x±,R), these two joint pdfs are related to each other, 



logio x 2 = log 10 xi + log 10 a . 



(11) 




(12) 



Pi R (x 1 ,R)=x 1 Pi2(x 1 ,x 2 ) 



(13) 
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where we use a modified ratio R = x 2 jx\ . From this relation, the detailed quasi-balance (12) 
is rewritten in terms of Pm(xi,R) as follows: 

Pir(x u R) =aR- 1 P lR ((^j 1 '\a 2 R- 1 ) . (14) 

Substituting the joint pdf Pir(x±, R) for the conditional probability Q(R\x±) defined in Eq. (4), 
the detailed quasi-balance is expressed as 

Pi(xi) aQ(a 2 R- 1 \(x 2 /a) 1 / 9 ) 

— (15) 



Pi((x 2 /a) 1/e ) R Q{R\xi) 

a QjgiR- 1 ) _ 
= R^Rj- = ° {a) • (16) 

We have used Gibrat's law in Eq. (15). 

By expanding Eq. (16) around R = a, the following differential equation is obtained 

a G'(a) 6 Pi (xi) + xi P'(xi) = . (17) 

The solution is given by 

Pi(xi) = Ci Xl - aG 'W 6 . (18) 

Here we consider two power-law distributions N±(> x\) cx xi^ 1 and A^(> x 2 ) cc x 2 ^ 2 , which 
lead 

Pi(xi) = Ci xf 1 " 1 , (19) 

P 2 (x 2 ) = C 2 xs"^- 1 . (20) 

From Eq. (18) and (20), we identify that a G'(a) 9 = ji\ + 1. On the other hand, by taking 
x\ — > {x 2 /a) 1 ' e transformation, Eq. (18) is rewritten as 

Identifying that C 2 = C\/a l l e and Pi((x 2 /a) 1 '' 9 ) = P 2 (x 2 ), we obtain a G'(a) = fj, 2 + 1. Conse- 
quently the relation between fii, fi 2 and 6 is expressed as 

^±1 = e . (22) 

This is the dynamical equation which non-equilibrium economic systems satisfy. We confirm 
that the empirical data satisfy this correlation in Fig. 12. 



4 Conclusion 

In this paper, we have investigated the dynamical behavior of non-equilibrium system in 
the high scale region by employing data on the assessed value of land in 1983 - 2005 Japan. 
By assuming the detailed quasi-balance (12) and Gibrat's law (3), we have derived a relation 
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between the change of Pareto index and a symmetry in the detailed quasi-balance. The relation 
(22) has been confirmed in the empirical data nicely. The result in this paper is considered to 
be applied to other long-term economic data. The equation (22) is the dynamical one which 
non-equilibrium systems satisfy. 

We should comment on two separations between 6 and (/xi + l)/(/X2 + 1) in Fig. 12. An 
abrupt jump of Pareto index between 1985 and 1986 (2001 and 2002) is observed in Fig. 6. This 
means that the system changes vigorously in this period, where the symmetry is not represented 
as Eq. (11). Nevertheless, the dynamical equation (22) is valid in almost all the other quasistatic 
periods. 

The finding in this paper is a first step to explain varying macro phenomena by using laws 
in non-equilibrium systems. For the next step, we should investigate the dynamical behavior in 
the middle scale region. 
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Figure 1: Cumulative number distribution iV(> x) of land prices in 1983 - 1985. The number 
of the data points is "6,262", "6,706" and "7,169", respectively. 
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Figure 2: Cumulative number distribution iV(> x) of land prices in 1986 - 1991. The number 
of the data points is "7,613", "8,120", 8,739", "9,321", "9,927" and "10,603", respectively. 
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Figure 3: Cumulative number distribution iV(> x) of land prices in 1992 - 1996. The number 
of the data points is "11,641", "15,538", "20,807", "24,753" and "25,091", respectively. 
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land price data in 1998 
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Figure 4: Cumulative number distribution iV(> x) of land prices in 1997 - 2001. The number 
of the data points is "25,845", "26,524", "27,073", "27,738" and "28,294", respectively. 
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land price data in 2003 
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Figure 5: Cumulative number distribution iV(> x) of land prices in 2002 - 2005. The number 
of the data points is "29,279", "30,138", "30,775" and "31,230", respectively. 
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Figure 6: Annual change of Pareto index ji from 1983 to 2005. 
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Figure 7: The scatter plot of all pieces of land assessed in the database, the values of which in 
1986 (x\) and 1987 fa) exceeded 10 4 yen/m 2 . The number of data points is "1,603". 
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LandPricex , (inyen/m 2 ) 

Figure 8: The scatter plot of all pieces of land assessed in the database, the values of which in 
1989 (an) and 1990 (x 2 ) exceeded 10 4 yen/m 2 . The number of data points is "3,886". 




10' 



10 5 10 6 10 7 

LandPricex , (inyen/m 2 ) 



10" 



Figure 9: The scatter plot of all pieces of land assessed in the database, the values of which in 
1983 (m) and 1984 (x 2 ) exceeded 10 4 yen/m 2 . The number of data points is "1,160". 
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LandPricex , (inyen/m 2 ) 

Figure 10: The scatter plot of all pieces of land assessed in the database, the values of which in 
2004 (xi) and 2005 (x 2 ) exceeded 10 4 yen/m 2 . The number of data points is "5,739". 
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Figure 11: Annual change of 6 and a of Eq. (11) in the year (x\,X2) = (1983,1984) 
-(2004,2005). 
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Figure 12: Annual change of 6 and (/ii + l)/(/i2 + 1) in the year (xi,X2) = (1983,1984) — 
-(2004,2005). 
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